Spin edge helices in a perpendicular magnetic field 
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We present an exact solution to the problem of the spin edge states in the presence of equal 
Bychkov-Rashba and Dresselhaus spin-orbit fields in a two-dimensional electron system, restricted 
by a hard-wall confining potential and exposed to a perpendicular magnetic field. We find that 
the spectrum of the spin edge states depends critically on the orientation of the sample edges with 
respect to the crystallographic axes. Such a strikingly different spectral behavior generates new 
modes of the persistent spin helix — spin edge helices with novel properties, which can be tuned by 
the applied electric and magnetic fields. 

PACS numbers: 72.25.-b, 72.15. Gd, 85.75.-d 



Introduction. Edge states are a defining factor in 
many prominent transport phenomena in condensed mat- 
ter physics. They emerge and are protected against scat- 
tering in the quantum Hall systems by applying a perpen- 
dicular magnetic field to a two-dimensional conductor, or 
are formed at the interfaces of topological insulators [l| 
by deforming the bulk band structure in the presence of 
strong spin-orbit coupling. The interplay of the three 
effects: the magnetic field quantization, the spin-orbit 
coupling, and the confinement by sample edges, is yet 
largely unexplored. 

Spin-orbit coupling is an important tool to manipu- 
late electron spins in solids by purely electric means 

In particular, the spin Hall effect or the spin 
Hall drag allow to create spin accumulation across 
the transport channels. The Bychkov-Rashba (BR) spin- 
orbit coupling , due to the structure inversion asymme- 
try, and the Dresselhaus (D) coupling |8| , due to the con- 
finement quantization of the bulk spin-orbit interaction 
(SOI) dominate in many semiconductor heterostructures. 
Both BR and D couplings can be tuned by electric gates, 
asymmetric doping, or strain, allowing for efficient spin 
control If the strengths of the BR and D couplings 

are equal, there can be long-lived, persistent spin helices 
(PSH) formed and protected against spin relaxation 13 1. 
This exciting phenomenon was recently observed 14l.ll5|. 

Here we study the complex interplay of the SOI, cy- 
clotron effects of an external magnetic field, and the 
hard wall confinement in a generic zinc-blende two- 
dimensional electron system (2DES) grown along [001]. 
We present an exact solution to the problem of the spin 
edge states for the equal BR and D SOI strengths. The 
spectrum of the edge states is strikingly different if the 
edges run along [110] or [110] directions. Depending on 
the relative sign of the BR and D couplings, in those 
two symmetry directions either (a) the spectrum is spin 
degenerate or (b) the spin splitting of the edge channels 
becomes maximal, i.e. either spin polarized channels and 
spin current oscillations jl3] are possible or not. We find 
that in (a) the non-Abelian gauge field via the built-in 



spin-dependent phase factor generates spin edge helices 
(SEH) with a precession angle that depends on the trans- 
verse distance from the edge. In (b) the shifting prop- 
erty of the spectrum allows the existence of SEH with a 
precession angle that depends on the distance along the 
propagation direction. In strong magnetic fields the pre- 
cession angle of the SEH is quantized in case (a), while 
a periodic helical structure, extended along the edge, is 
produced in case (b). In weak magnetic fields we find 
interesting new spin resonances when the cyclotron mo- 
tion is commensurate with the spin precession [l7| . Ex- 
perimentally, a strong reduction of spin scattering rate 
towards the sample edges should be observed. 

As an important application of our theory we propose 
extending the experimental setup in Ref. [ij by exposing 
the 2DES additionally to a perpendicular magnetic field 
B. According to the experimental findings from Ref. Il4. 
the enhancement of the PSH by about two-order of mag- 
nitude occurs only at intermediate temperatures about 
T ^ 100 K. Its rapid drop with lowering of T shows that 
the spin Coulomb drag |18l420l | via a strong increase of 
the diffusion coefficient destroys the PSH enhancement, 
for which a finite momentum Q is needed in _B = 0. In 
the case of finite B we find that in order to excite SEH 
in the (b) configuration a finite momentum shift is still 
needed between the spin components but in the (a) con- 
figuration, where the energy is degenerate, the SEH ex- 
ists for Q = 0, i.e. without momentum difference along 
its propagation direction. This should reduce the role 
of spin Coulomb drag in suppressing the enhancement 
of SEH, which can be useful for spintronic applications. 
Notice that in three-dimension a nonequivalence in (110) 
and (110) has been already observed by means of the 
electron paramagnetic resonance and the electron-dipole 



17, 21 



spin resonance 

Theoretical concept. We consider electrons in a 2DES 
exposed to a perpendicular magnetic field B along [001]. 
The electrons are additionally confined by an infinite po- 
tential V (x) = oo for a; < 0. Then in the presence of the 
BR and D SOI the Hamihonian is H = Ho + Hsoi + V{x) 



2 



1-5 «: 



[1 10] 




FIG. 1. (Color online) Two physically difTerent configurations 
of the 2DES, restricted with a hard- wall confining potential 
at a; = (filled areas). The spin edge states propagate in 
the skipping orbits along [110] and [110] axes, respectively, 
parallel and perpendicular to the direction of electron spins. 
The Fermi contours in the absence of the magnetic field are 
shown in the momentum plane with arrows indicating the 
directions of spin. 



where Hq — 7r'^/2m* desribes a free particle in a quan- 
tizing magnetic field, m* denotes the electron effective 
mass and n = p — (e/c)A the kinetic momentum with 
p — —ihV; the electron charge is — e. There are two 
preferential directions in the plane of 2DES: i) the 

direction of the sample boundary along which edge states 
propagate, and ii) the direction of the electron spin in the 
presence of BR and D SOI of equal strength, determined 
by the crystallographic axes. The relative orientation of 
these two directions determines two distinct configura- 
tions, shown in Fig. [TJ In these configurations we choose 
the (x, y) coordinate system such that the sample bound- 
ary is always along y. Then in the configuration (a) in 
Fig. [U in which the axes x and y are along [110] and 
[TlO], we have 



HsOI = ("i? " cud) T^y^x - {OLR + ao) TTx^y 



(1) 



In the configuration (b), the axes x and y are along [TlO] 
and [TTO] and the Hamiltonian is 



HsOI = ("i? + Oio) TTy&x - {ur - ao) T^x^y 



(2) 



Here an and an are, respectively, BR and D spin-orbit 
coupling constants and crj.,cTy are the Pauli matrices. 
We consider magnetic fields which are strong enough to 
quantize the electron spectrum but weak enough to cause 
much smaller Zeeman splitting than the SOI induced en- 
ergy splitting 1^ 2^ . 

We choose the Landau gauge with A{x) = (0,.Ti?, 0) 
and making the ansatz 



— pikyy 



Xky{x) , f = (a;,?/) , 



(3) 



reduce the two-dimensional Schrodinger equation to the 
effective one-dimensional problem in the x direction. The 



transformed Hamiltonian U^HU, generated by the global 
unitary transformations 



jja ^ 



V2VI 1 



and U" = 



V2 \ I I 



(4) 



becomes diagonal in the (a) and (b) configurations so 
that in the case of a/? = ajj—a the wave functions Xky (x) 
satisfy the following equation 



- — h i^aa- 
ax 



^+2 



1 {x- X (ky) - "fba^Y 



Xky (x) = , 
(5) 



Here the coefficients a and b are given by 



a = 1, b = 0, for the (a) configuration , (6) 
a = 0, b = 2, for the (b) configuration . (7) 

In Eq. = i' + -f'^ and we express the total electron en- 
ergy E ^ {ly + 1/2) hojB and the length x — )• xIb/V^ in 
units of the cyclotron energy, Hujb — heB/m*c, and the 
magnetic length, Ib = \Jhc/eB. The dimensionless SOI 
coupling constant 7 = ^/2a/vB, where the cyclotron ve- 
locity vb = h/m*lB, while the dimensionless coordinate 
of the center of orbital rotation X (ky) = y/2kylB- 

In the (a) configuration the SOI with equal BR and 
D strengths induces a finite non-Abelian gauge field that 
depends on the spin orientation, keeping the cyclotron 
rotation center unshifted. In contrast, in the (b) configu- 
ration the SOI induces only a spin-dependent shift of the 
cyclotron rotation center. 

With the unitary transformation 



Xky (x) = cxp (-170x6-2) (l)k^ (x) 



(8) 



we can eliminate the non-Abelian gauge field and map 
the SOI problem of the equal BR and D strengths to the 
SOI free, shifted edge state problem as follows: 



(x - S76) (f>i (x) = 



Here we introduce the operator 



hiy (x) 




{x-x{ky)y 



(9) 



(10) 



and s = ±1 labels the spin f and ^ states in the new spin 
basis, created by The general solution of Eq. [5] is 
given in terms of the parabolic cylindric functions I?^ (x) 
so that the spin edge states are given by 

ipl (r) — exp {—isjax + ikyy) [x — X [ky) — sjb) . 

(11) 

For sufficiently large positive X (ky) the solution pT|) 
recovers the exact spectrum of the dispersionless bulk 
Landau levels, E^d^) = + ^ — 7^) Swb with the index 
I = 0, 1, 2 . . . [23j; this is valid in both configurations since 
the index p, is the same. The shift of all bulk Landau lev- 
els due to SOI is independent of spin so they remain spin 
degenerate. 
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FIG. 2. The energy spectrum of spin edge states in the 
presence of Bychkov-Rashba and Dresselhaus SOI of equal 
strengths. The dashed and solid curves correspond to the (a) 
and (b) configurations in Fig. 1. 



Energy spectrum. The energy spectrum of the 
spin edge states is obtained by letting the wave 
functions vanish at the sample boundary: 

D^(^^,y-^{x~X{ky)~sjb)\^_^ — 0. The calculated 
energy branches, Esi{ky,^), are shown in Fig. [5] As seen 
in the (a) configuration the energy spectrum remains 
spin degenerate for all values of ky so that the spin 
edge states at the Fermi energy, Ep, are not separated 
in space. The wave functions (jlip contain finite spin- 
dependent phase factors (a = 1), which describe the 
precession of spins with the opposite rotations along the 
electron propagation axis y. 

In contrast, the energy spectrum in (b) develops, for 
a given principal quantum number two spin branches. 
The degeneracy is lifted already for positive X {ky) ^ 2; 
the spin-splitting increases strongly with \X [ky) \ for 
negative X{ky). The spin-dependent phase factor in 
([TT|) vanishes identically (a = 0) and the effect of equal 
BR and D SOI on the spin edge states is reduced to the 
spin-dependent shift of the energy branches in the mo- 
mentum space. At certain values of Ep there exist two 
spatially separated current carrying states with opposite 
spins, whose directions are locked globally by the geom- 
etry of the configuration (b). 

Thus, the spectrum in Fig. [2] shows that the exter- 
nal magnetic field creates edge states along the [110] and 
[110]. The linear motion of the edge states via the BR 
and D SOI of equal strengths induces an effective mag- 
netic field along [110] and [110], respectively, that is ei- 
ther perpendicular or parallel to the direction of the spins 
thereby keeping the spin degeneracy or resulting in the 
maximal spin-splitting, respectively, in (a) or (b). 

Spin edge helices. The energy bands of the spin edge 
states possess a shifting property along y. 



where the shifting wavenumber Q = y/2'yb/lB is finite in 
(b) and zero in (a) . As in the zero field [13| , we introduce 
the following operators in the transformed spin basis 

Sq (rl = < (r)^ <.Q (rO , 5+ (r) = {r)^ (r) 

(13) 

^0 (^1 = V-I, i^^ < - < (r)^ < (0 , 

which commute with the system Hamiltonian owing to 
the shifting property (|T2]) and the SU(2) symmetry of 
the system Il3| . The non-diagonal operators 



exp {zL2iajx ± ibjy) [x ~ X{ky) + 67] 



Eii{ky,-f) = E^i {ky 



,7), 



(12) 



(14) 

represent long-lived .spin edge helices of 2DES in the pres- 
ence of a perpendicular magnetic field. Here the coordi- 
nate y is also dimensionless, y — )■ ylB/V^- As seen from 
Eq. [T3J in contrast to the spin edge states, the helical 
edge modes have finite spin- dependent phase factors in 
both configurations. Going back to the initial spin basis 
we see that in the (a) configuration (a = 1 and b — 0) 
the factor exp {±2i^x), existing also for the edge states, 
describes the spin precession in {x, z) plane with the pre- 
cession angle d{x) = A'yx, depending on the transverse x 
coordinate, along which the helices are confined by the 
magnetic field via the factor Dfi [x — X{ky) -\- 67]^. On 
the contrary, in the (b) configuration {a — Q and b — 2) 
the factor exp (±2i7y) is inherent only to the SEH. This 
factor arises from the plane wave functions in Eq.[3]due to 
the shifting property (|12p and describes the spin preces- 
sion in (y, z) plane with the precession angle 'd{x) = 47?/, 
depending on the y coordinate along the free propagation 
direction of the SEH. Thus the combined effect of the 
perpendicular magnetic field and the confining potential 
on the SEH critically depends on the orientation of the 
edges of 2DES relative to the crystallographic axes; as we 
will see it depends also on the strength of the magnetic 
field. 

In quantizing magnetic fields the spin helices are 
strongly localized along the transverse direction x around 
its average position xi {ky). In the limit of large nega- 
tive ky, xi {ky) is approximately independent of fcj, and 
takes discrete values xi in the different channels I |16l.l27l|. 
Therefore, in the (a) configuration the precession angle 
is quantized around di = 'i'yxi for the spatially separated 
edge channels I (see bottom of Fig. 3. a). Meanwhile in 
the (b) configuration the precession of spin by an angle 
•ffi = 47^ develops a spatially periodic structure along 
the y direction (see bottom of Fig. 3.b), similarly to the 
usual PSH in an infinite 2DES in the magnetic field free 
case [13, 24-2^. Thus, by switching between the (a) and 
(b) configurations, i.e. by tuning the BR and D coupling 
strengths either to the an = an or an = —ao case, one 
can realize a selection mechanism between the two alter- 
natives of the quantized spin edge helices and of the free 
spin edge helices. 
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FIG. 3. Spin edge helices in a perpendicular magnetic field. 
Figs, (a) and (c) (Figs, (b) and (d)) correspond to the (a) con- 
figuration (the (b) configuration). Figs, (a) and (b) (Figs, (c) 
and (d)) correspond to the limit of strong (weak) magnetic 
fields. The bottom (top) of each figure represents the heli- 
cal structures in the edge channels (in the quasibulk Landau 
states). 



In the limit of large positive fcj,, the electrons are lo- 
calized in quasibulk Landau states around their rota- 
tion center X(ky) ^ 1, which increases linearly with 



ky [H, 27 1 . However, for a given ky one can neglect 



weak oscillations of the spin precession angle since 
in strong fields Ib is much smaller than the spin-orbit 
length, Aso ~ l/2ma, and the variation of "di during 
the period of the cyclotron rotation is of the order of 
7 = 2malB ^ 1- Thus, the precession angle in the con- 
figuration (a), di « A'yXiky), is independent of y and 
varies only along x with ky (see top of Fig. 3a) while in 
the (b) configuration, w 477/ is independent of x but 
varies along y (see top of Fig. 3b) . 

In weak fields and for Z 1, the quasi-classical de- 
scription is valid and the x coordinate oscillates in time as 
X = X {ky) + ^/2v + Icosujut. For the edge states X iky) 
can be negative, —00 < Xg {ky) ^ ^/2v + 1 so that in the 
(a) configuration the spin precesses along x by an an- 
gle ^ ^ 2do where i?o = '^l\/2v + 1. Meanwhile for 
the bulk Landau states X [ky) ^ ^J2v -\- 1 and the preces- 
sion angle in the (a) configuration varies within the range 
i?s — i?o ^ ^ ds + Therefore, in the magnetic fields 
and electron concentrations such that A'yy/2vF + 1 = 7r/2 
with 2vp + 1 = Ep/huB, the spin resonance effect takes 
place as shown in the top of Fig. 3c: the spin makes 
one full precession each time when the electron makes a 
full cyclotron rotation. Recently, the spin resonance phe- 
nomenon, driven by the SOI and perpendicular magnetic 
fields, has been experimentally demonstrated in Ref. [2^ 
by observing a strong suppression of spin relaxation. As 
seen in the top of Fig. 3d and 3c the spin resonance for 
the Landau states in the (b) configuration differs from 
that in the (a) by a simple rotation of the cyclotron orbit 
because in the (a) = Ajy depends only on y while in 
the (b) t} = 4'yx depends only on x. 

For large negative ky the physical picture of the spin 
resonance in the (a) and (b) configurations differs essen- 
tially (c/. the bottom of Figs. 3c and 3d). In (a) = 47a; 



depends only on the x, along which the electron motion 
oscillates within a finite range. Therefore, the preces- 
sion angle is independent of the electron motion along 
the edge. In contrast, in (b) ^ = Ajy depends on the 
y coordinate, along which the electron propagates freely. 
Therefore, in (b) the spin precession generates a spatially 
periodic structure of the SEH along y and as seen in the 
bottom of Fig. 3d in the resonance case the spin changes 
its sign each time the electron makes a half circle in its 
skipping orbit. 

In conclusion, we present a theory of persistent spin 
edge helices, which exhibit novel features, tuneable by 
electric and magnetic fields. We show that either a peri- 
odic structure of spin edge helices along the sample edges 
or a helical structure with a quantized precession angle 
along the transverse direction is realized. 
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